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ABSTRACT 

The next generation of cosmic microwave background (CMB) experiments, galaxy 
surveys, and high-redshift observations can potentially determine the nature of the dark 
matter observationally. With this in mind, we introduce a phenomenological model 
for a generalized dark matter (GDM) component and discuss its effect on large-scale 
structure and CMB anisotropies. Specifying the gravitational influence of the otherwise 
non-interacting GDM requires not merely a model for its equation of state but one for 
its full stress tensor. From consideration of symmetries, conservation laws, and gauge 
invariance, we construct a simple but powerful 3 component parameterization of these 
stresses that exposes the new phenomena produced by GDM. Limiting cases include: a 
particle component (e.g. WIMPS, radiation or massive neutrinos), a cosmological con- 
stant, and a scalar field component. Intermediate cases illustrate how the clustering 
properties of the dark matter can be specified independently of its equation of state. 
This freedom allows one to alter the amplitude and features in the matter power spec- 
trum relative to those of the CMB anisotropies while leaving the background cosmology 
fixed. Conversely, observational constraints on such phenomena can help determine the 
nature of the dark matter. 

Subject headings: cosmology: theory - dark matter - large-scale structure of the universe 
- cosmic microwave background 
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1. Introduction 

Upcoming cosmic microwave background (CMB) mis- 
sions, galaxy redshift surveys, and high-redshift obser- 
vations will produce such a wealth of high-quality data 
that even the extended cold dark matter (CD M) model 
with 11 free parameters (e.g. Jungman et al. 1996) may 
fail to fit them. One must face the very real possibility 
that none of our current ab initio models will survive 
the upcoming confrontation with the data. 

How might one generalize the CDM model? The cor- 
nerstone of all modern cosmologies is, and will likely 
remain, gravitational instability in a world-model that 
is homogeneo us and isotropic on the large scale (e.g. 
Peebles et al. 1991). Models for the dark matter sector, 
on the other hand, are presently limited by the number 
of candidates that are considered well-motivated from 
the particle physics standpoint, e. g. WIMPS (CDM), 
massive neutrinos, scalar fields (see |Coblc, Dodclson"fc" 
Frieman| 19 96; Ferreira fc Joyce- 1997; paldwell, Dave fc 
stcinhardt 1998, for recent assessments) and topological 
defects QSpergel fc Perj 1997). If none of these candi- 
dates survive the confrontation with high-precision cos- 
mological measurements, we will be forced to solve the 
inverse problem: can one determine the nature of the 
dark matter and reconstruct the model for structure 
formation directly from the observations? 

Explorations of the dark matter sector have been un- 



dertaken recently by Turner & White (1997) and Cald 
well, Dave & Steinhardt (1998) who considered dark 
matter with an arbitrary equation of state that pos- 
sesses no fluctuations and scalar-field type fluctuations 
respectively. The former case does not present a com- 
plete theory of structure formation as it can only apply 
below the horizon at any given time for adiabatic mod- 
els. The latter case provides an interesting example of 
an exotic dark matter component but does not exhaust 
the possibilities for its gravitational properties. For ex- 
ample, a hot dark matter component has an equation 
of state like matter today but clusters only up to a fi- 
nite scale. This scale is substantially below the current 
horizon scale, where it would be for the analogous scalar 
field model ( Ferreira fc Joyce] 1997). 

Starting from the general principles of symmetry, 
energy-momentum conservation and gauge invariancc, 
we build in ^ a phenomenological parameterization of 
the dark matter that includes all of the currently popu- 
lar dark matter theories as special cases. The main re- 
sult is the link established between the clustering prop- 
erties of the dark matter and the model for its underly- 



ing stresses. Importantly, these properties are not de- 
termined by the background equation of state in the 
general case. Through their effect on the growth rate 
of perturbations, discussed in §|], these clustering prop- 
erties manifest themselves themselves as independent 
features in the power spectra of large-scale structure 
and CMB anisotropies as discussed in ^ and §||. We 
summarize the main results in §[| and present a short 
Appendix that highlights the scalar field case. 

2. Dark Matter Properties 
2.1. General Principles 

The gravitational influence of an arbitrary dark mat- 
ter component is controlled by its stress-energy tensor 
r Ml/ (x, rf), where 77 = J dt/a is the conformal time with 
a(t) as the scale factor normalized to unity today. In 
general, the symmetric 4-tensor T^ v has 10 components 
which can be divided into 4 classes: the energy den- 
sity p g (1), the isotropic stress or pressure p g (1), the 
momentum density (p g + p g )v l g (3), and the anisotropic 
stress p g 7T l J (5). The 5 components of the anisotropic 
stress can be further separated by their transforma- 
tion properties under rotations into a scalar component 
(1), vector components (2), and tensor components (2). 
Since only scalar components exhibit gravitational in- 
stability, we hereafter neglect the vector and tensor con- 
tributions]^ Energy-momentum conservation T^ u ' u = 
introduces 4 constraints, leaving only two independent 
parameters for the dark matter. One can choose these 
to be the press ure p Q an d scalar anisotropic stress am- 
plitude Ti g [see Bardeen 1980, Eq. (2.18)] without loss 
of generality. 

The isotropy of the background implies that the 
anisotropic stress and momentum density can only be 
present as a perturbation. The conservation laws then 
reduce to a single relation, 



Pg 
Pg 



-3(1 + w g )- 



(1) 



where overdots represent conformal time derivatives 



and 



G^ u = reduce to 



p g /p g - Likewise the Einstein equations 



iirG 9 \ -s 
3 ^ 



> pi — K , 



(2) 



2 Vector and tensor contributions do however affect CMB 
anisotropies and can act as additional degrees of freedom when 
normalizing large-scale structure to CMB measurements (e.g. 
|Caldwcll fc Steinhardt] 1998). 
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where the sum is over the density contributions of all 
matter species and the background curvature is K — 
—Hq(1 — fitot) with fitot = Yli^i- As usua l the ex- 
pansion rate today is given by the Hubble constant 
Hq = {a/a) a= i — 100ft- km s _1 Mpc -1 , to which each 
species contributes according to its fraction of the crit- 
ical density fij = 8nGpi/3HQ. 

The background evolution is thus completely speci- 
fied by the equation of state w g (a). This is not true of 
the perturbations, where we are left with the freedom 
to specify Sp g and 7r 9 . It is convenient to separate out 
the non-adiabatic stress or entropy contribution 



PgTg = S Pg ~ CgS P g , 

where the adiabatic sound speed is 



Pg_ 
Pa 



1 



3 1 



W 



(3) 



(4) 



Therefore, p„ = 



' J gPa 



does not imply Sp g = w g Sp g (due 



to temporal or spatial variations in w g ) and further- 
more if V g ^ 0, the function w g (a) does not completely 
specify the pressure fluctuation. 

In the following, we adopt the notation of Hu et al. 
(1998). For brevity, we present only the aspects of per- 
turbation theory that are altered by the presence of 
GDM. Energy-momentum conservation yields the con- 
tinuity equation for the density fluctuation S g = Sp g / p g 



(kvg + 3h 5 )-3--^Tg, (5) 



1+Wg 

and the Euler equation 



v g = -«(l-3c>)v g + —M-k6 g + —JL-kr a 



a W 9 • ! • w g ' l + wg 

2 w„ 



-(1 -3K/k 2 )kir g + kh v . 



(6) 



3 1 + Wg 

The metric sources h$ and h v depend on the choice of 
gauge and are 



h s = 



h v 



$ 




Synchronous, 
Newtonian, 

Synchronous, 
Newtonian, 



(7) 



Note that = h/6 in the notation of Ma fc Bertschingci 
(1995). Seed perturbations (e.g. defects) whose total 
contribution is first order in the perturbations can also 
be modeled in this manner by rewriting equations (|l]), 
(||) and (^|) in terms of Sp g and (p q + p g )v instead of 



the relative perturbations (see eq. |A9| ) 



2.2. Stress Model 

Up to this point, we have made no assumptions what- 
soever about the nature of the dark matter since w g , T g 
and TT g have been left as free functions. We now need to 
parameterize these functions. Dark matter with w g < 
is favored by current observational constraints, such as 
the combination of the ages of globula r clusters and the 
high Hubble constant measurements (Chiba, Sugiyama 
|fe Nakamura| 1997; |Caldwcll, Dave fc Steinhardt| 1998) 



as well as supernova luminosity distance measures (Perl- 
mutter, et al. 1998;Riess, et al. 1998). If w g < and is 
slowly- varying compared with the expansion rate (a/a) 
such that c 2 < 0, the adiabatic pressure fluctuation 
produces accelerated collapse rather than support for 
the density perturbation. In a GDM-dominated uni- 
verse, perturbations would rapidly go nonlinear once 
the sound horizon has been crossed \k J c g dn\ > 1. This 
situation is unacceptable for a model of structure for- 
mation. 

In this w g < regime, it is interesting to con- 
sider whether non-adiabatic pressure can act to stabi- 
lize the perturbation. This requires a relation of the 
type WgTg oc 6 g . One is not however allowed complete 
freedom in establishing this relation. Adiabatic pres- 
sure, density, and velocity perturbations are gauge de- 
pendent whereas non-adiabatic pressure perturbations 
are not. Therefore, stabilization in one frame of refer- 
ence does not equate to stabilization in another. One 
should avoid having the properties of the GDM depend 
on an arbitrary choice of frame and hence unphysically 
on the perturbations in the other species. This require- 
ment can be achieved by defining the relation in the rest 
frame of the GDM where T°, = 0, 



v 'a 



(8) 



-gig 

We further assume the effective sound speed c^ s is only 
a function of time. If c^ ff > 0, pressure support is 
obtained. The gauge transformation into an arbitrary 
frame gives 

Of 3 ^ = Sg + 3^(1 + Wg)(Vg - B) / k , (9) 

yielding a manifestly gauge-invariant form for the non- 



adiabatic stress (see Bardeen 1980; Kodama & Sasaki 
1984). Here B represents the time-space component 
of metric fluctuations and vanishes in both the syn- 
chronous and Newtonian gauges. The Euler equation 
can then be rewritten as 



a 

-~V n 



biff^rest) 



W 



3 



-kh - -AKik 2 y g + hh v (in) 



3T - 

Thus Cg ff may be thought of as a rest-frame sound speed. 
By inspection of equations ( p^| ) and (|5|) , we determine 
that the critical scale for stabilization is the effective 
sound horizon, 



Spergel & White 1997 for a parallel treatment). Note 



Soff 



di]c cS . 



(11) 



This assumes w g is not varying at a rate much greater 
than the expansion rate (see Appendix). 

The anisotropic stress can also affect the density 
perturbations. A familiar example is that of fluid, 
where it represents viscosity and damps density per- 
turbations. More generally, the anisotropic stress com- 
ponent is the amplitude of a 3-tensor that is linear in 
the perturbation. A natural choice for its source is kv g , 
the amplitude of the velocity shear tensor d l v g . How- 
ever it must also be gauge invariant and generated by 
the corresponding shear term in the metric fluctuation 
Ht- The relationship between velocity/metric shear 
and anisotropic stress may be parameterized with a 
"viscosity parameter" c 2 , 



- 2 . 

-"vis 



3-TTg ) = 4c^ is (fcWg - H T ) 



where in the Newtonian gauge Ht = and in the 
sy nchronous gauge Ht — hr{= —h/2 — 3-q, Ma fc 
Bei tschingci 1995). The specific form of this equation 
is designed to recover the free streaming equations of 
motion for radiation with an approximate closing of the 



angular moment hierarchy at the quadrupole (Hu et al 



1995). The physical interpretation of equation ([12]) is 
that the anisotropic stress will act to damp out velocity 
fluctuations on shear-free frames (Ht — 0) if c^ is > 0. 
We call s v i s = J Cyisdrj the viscous scale. 

We shall see in the following sections that this pa- 
rameterization captures many of the essential features 
of GDM as it corresponds to a means of altering its 
clustering properties. In the limit that (w g , Cg ff , c^ is ) — > 
(w g , 1,0) scalar-field dark matter is recovered exactly 
(see Appendix). CDM is similarly recovered with (0, 0, 0), 
radiation can be modeled as (1/3, 1/3, 1/3) and hot dark 
matter (HDM) or warm dark matter in a similar matter 
described in § )4.3) 

Cases that this ansatz does not cover involve mainly 
models in which stress fluctuations are not derived from 
density and velocity perturbations and so act as exter- 
nal sources for the perturbations. This situation occurs 



in the case of non-linear seed perturbations (see Hu 



that in such models both vector and tensor stresses must 
also be modeled to yield a complete theory for structure 
formation (see [Hu fc White) 1997b; )Turok, Pen, fc Seljak 
1997). 
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Fig. 1. — Large scale perturbation evolution (ks e g 3> 1) 



with GDM of 



-1/6, 



- cir 



1 (solid, scalar fields) 



and (dashed, stress-gradient free). The Newtonian 
(12) curvature $ is independent of c^ s and varies only when 
the background equation of state changes at a eq and 
the ratio of density perturbations in the 
GDM and matter depends on Cg ff . The cosmological 



a g . However, 



parameters here are fitot 
h = 0.7. 



i, n g = 0.9, fW = 0.0125, 



Equations (§), (§), and @ can now be introduced 
into any of the standard codes that solve the Einstein- 
Boltzmann equations. We employ the code of Hu et al 
(1995) for the examples below. We furthermore show 
adiabatic models for illustrative purposes but of course 
isocurvature models can be similarly obtained through 
a change in the initial conditions. 

3. Perturbation Growth 

By introducing a means by which fluctuations in the 
GDM are stabilized, we change the growth rate of fluc- 
tuations in the baryons and any CDM that may be 
present. The clustering properties of the GDM thus 
have consequences for both large-scale structure and 
CMB anisotropics as we shall see in §)i| and ||. Here we 
summariz e results for the g rowth rate of perturbations 
proven m |Hu fc EiscnstcirJ (1998). 
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3.1. Above the Sound Horizon 

Above the stabilization scale of all species and in the 
absence of background curvature (K — 0), perturba- 
tions evolve so as to keep the Newtonian curvature $ 
and potential W constant 

i 

fe 2 (# + $) = -87TGa 2 J2 Pi ir t , (13) 

except during periods when the dominant equation of 
state w change s. Here the potentials vary mildly such 
that ( |Bardccn| 1980) 



2 1 

<f> * = const. 

31 + w 



(14) 



For reference, dominance of GDM relative to matter 
occurs at 



l/3ii 



(15) 



The result is that during periods when the equation 
of state is slowly- varying, the total density fluctuation 
§(test) _ Pi/Y,Pi grows as 



5 {r 



1+3™ 



(16) 



We display an example of this large-scale evolution in 
Fig. | 

Since CMB anisotropics are only dependent on the 
time evolution of the potentials, the effect of GDM here 
is very weak. The exception is the w g — ► — 1 limit, where 
the potentials decay to zero and the cosmological con- 
stant case is recovered. Note however that a curvature 
term in the background and a dark matter component 
with w g — — 1/3 give identical contributions to the ex- 
pansion rate but are not similar in their contribution to 
large-angle CMB anisotropics. 

Finally, the division of density fluctuations between 
the GDM and the matter components, i.e. S Q /S m , de- 



per ds on the form of the GDM stress. This is because 



of the appearance of the non-adiabatic stress term in 
the continuity equation (^]) which relates the density 
and metric fluctuations. As c 2 s increases, S g /S m de- 
creases. This affects the amplitude of the matter power 



spectrum as we shall see (§4.1) 



10 2 



10 1 



xs 10° 



Gio- 1 



io- 



- 1 1 II II II 1 1 1 1 1 1 1 1 1 1 1 


Mill 


\ 1 1 II ± 










, ' '' ' 


= $x 10 






E A (rest) /T 
























^k=10- 2 Mpc- 1 

' 1 1 




T>v, mi 


io- 3 io- 2 


io- 1 




a 







io- 



Fig. 2. — Small-scale perturbation evolution for the 
GDM and the matter density perturbations in the same 
models as Fig. [l] with c 2 s = 1 (dashed lines) and 
Cg ff = (solid lines) . GDM perturbations stabilize once 
fcs c ff > 7r and their relative absence (Sg cst ^/(5m CSt ^ <C 1) 
then slow the growth of matter perturbations once the 
expansion is also GDM dominated a > a g leaving the 
potential $ to decay. 

3.2. Below the Sound Horizon 

Below the effective sound horizon of the GDM, its 
perturbations stabilize. If GDM dominates the energy 
density, Sg cst ' ) oscillates with a decaying amplitude of 



A oc a(- 1+3 ^)/ 2 c- 1/2 



(17) 



and rapidly becomes a smooth density component com- 
pared with the fluctuations in the other species. We 
display a case where there is also a non-negligible CDM 
component in Fig. 0. Once GDM dominates the en- 
ergy density of the universe at a g , then the smooth- 
ing of the GDM-component will also slow or halt the 
growth in the matter species. In particular, the growth 
will slow to a halt if w g < and scale as a p with 
Ap = v/l + 24f2 m /f2tot — 1 if w g = 0. Analytic so- 



lutions for how this process occurs are given in Hu & 



Eiscnsteinl (1998). 

The net effect is that if there is enough CDM to 
ensure that the universe was matter dominated some- 
time in the past, perturbation growth is suppressed 
by a scale-independent factor below the effective sound 
horizon s e ff at GDM-domination. This suppression de- 
creases until it disappears above s c s today. For the 
CMB, the suppression of growth in the density pertur- 
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bations causes the potentials to decay to zero, leading 
to a contribution potentially 36 times larger than the 
Sachs- Wolfe effect in the anisotropy power spectrum. 
We shall see in §5.2 why this limiting value is never 



reached in practice. 

On the other hand, if matter never dominated the 
expansion as is the case if CDM is absent, even more 
dramatic effects occur. In this case, even for c 2 s = 

0, growth is highly suppressed on small scales due to 
an extended period of radiation domination (assuming 
w g < 0). The controlling scale is therefore the horizon 
at GDM-radiation equality 

1. e. 



v 3eq 



—n g H (^-y 



(18) 



(19) 



Above this scale, perturbations grow as in equation (|16|); 
below this scale, perturbations only experience signifi- 
cant growth after GDM-domination. Because GDM- 
domination for w g < is delayed compared with an 
equivalent w g — universe, the critical scale is larger in 
such single component GDM models than in the CDM 
case. 

3.3. Viscosity 

Finally the viscous stresses of equation ( |l2| ) can dis- 
sipate fluctuations in the GDM. In Fig. ||, we show an 
example with w g — 1/3 GDM replacing the three mass- 
less neutrinos in an otherwise standard CDM universe 
(sCDM, O tot = 1 ~ On, h = 0.5, n b h 2 = 0.0125). 
The real neutrino background radiation contains an 
anisotropic stress due to the quadrupole moment of 
its temperature distribution. Modeling the neutrinos 
as GDM allows us to explore the consequences of its 
anisotropy by varying c^ is and also illustrates the phe- 
nomenological manifestations of the viscosity parame- 
ter. 

In the absence of anisotropic stresses (dashed lines, 
c vis — 0), perturbations in the GDM oscillate. Chang- 
ing c^ is to 1/3 to approximate the radiative viscosity of 
the real neutrinos, the density perturbations damp once 
fcsvis J> 7T- Note however, its effect on the gravitational 
potentials $ and well after sound horizon crossing 
when large-scale structure is observed is negligible. This 
is because the pressure fluctuations are sufficiently ef- 
fective to make the GDM perturbations smooth in com- 
parison to the growing species (see Fig. ^). The extra 




Fig. 3. — Viscous effects with GDM w g = 1/3 replacing 
the three species of massless neutrinos in the sCDM 
model (Qtot = 1 ~ O m , h = 0.5, n b h 2 = 0.0125). With 
the viscosity parameter set to mimic radiation c^ is = 
1/3 (solid lines) the perturbations in the GDM decay 
whereas with c^ is = 0, they do not. This distinction 
has a negligible effect on the behavior of the potentials 
$ and 'J well after sound horizon crossing 



smoothing due to viscous damping affects perturbations 
little. 

The anisotropic stress does change the behavior of 
the potentials at early times because it enters directly 



in to the Poisson equations (O). We shall see in §5.3 



that correspondingly viscous effects are more important 
for the CMB than large-scale structure so long as \c 2 s | > 

I _2 I 



4. Large-Scale Structure 

The large-scale structure of the universe depends on 
the detailed properties of the GDM. The main result 
is that the clustering scale becomes independent of the 
equation of state of the dark matter. By changing the 
growth rate of perturbations, the clustering properties 
change the amplitude of and features in the matter 
power spectrum. Here we present concrete examples 
of this process that include scalar fields, radiation, and 
hot dark matter as special cases. 

4.1. Sound Horizon and Scalar Fields 

The introduction of a stabilization scale or effective 
sound horizon s c g for the GDM places a feature in the 
matter power spectrum between that scale at GDM- 
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domination and today, i.e. 



s eff( a = !)<*< s ff( a = a g)- 



(20) 



The limiting cases are the scalar-field example, where 
the effective sound speed is the speed of light 



= 1 



= 0. In 



cir 



and the pressure-gradient free case with qf ff 
Fig. H we display the effect of varying c 2 s on the power 
spectrum of the matter holding the background cosmol- 
ogy fixed. These models have been consistently nor- 
malized to the COBE CMB anisotropy measurement at 
large scales via the fitting form of Bunn fc White] (1997) 
(their equations [17]-[20]). 




10- 2 

k (h Mpc~ 



Fig. 4. — The effective sound horizon and the COBE 
normalized matter power spectrum. Raising c 2 s from 
to 1 (solid lines) introduces a feature between the 
effective sound horizon at GDM-domination and that 
scale today. Here a s — (0.75,0.29,0.25) and o^o/ag = 
(0.16,0.17,0.16) for c 2 cS = (0,1/6,1). These models 



have w a 



-1/6, 



0, n t 



1, n n 



0.65, 



flbh 2 — 0.0125, h = 0.7 and tilt n = 1. For comparison, 



anisotropics such that the COBE signal drops relative 
to the potential fluctuations. These two effects com- 
pete such that the change in normalization of P(k) at 
large scales with c 2 s is non-monotonic. The shape of the 
power spectrum, above and below the transition region 
of equation ( p0| ) remains that of a CDM model with the 
same f2 m /i and Q h 2 (c.f. (t§o/(t% in Fig. Q). 



^(a) Radiation 




0.001 0.01 0.1 1 

k (h Mpc- 1 ) 

Fig. 5. — (a) Modeling radiation. Shown here is the 
power spectrum for the model of Fig. || where CDM 
of uu g = 1/3 replaces the neutrinos of sCDM. Altering 
the viscosity parameter from c 2 is — 1/3 to has little 
effect on the power spectrum although 1/3 is a some- 
what better approximation at large scales, (b) Modeling 
HDM. The features of the mixed dark matter are well- 
reproduced by CDM with the same equation of state 
and c 2 ia = w a . The parameters here are sCDM with 



the corresponding A-model (w g -> -1, same parame- ^ = Q 2 placing part of the CDM 



ters, erg = 1.1, (T50/0-8 = 0.16, dashed lines), which fits 
the current large scale structure data, is also shown. 

The amplitude of fluctuations at the fiducial 8/i~ 1 Mpc 
scale, eg, is affected in three ways. The presence of a 
clustering scale reduces the growth rate below it leading 
to a relative suppression of small-scale power. However 
the absolute amplitude of large-scale fluctuations also 
changes with the clustering scale. As shown in the last 
section, decreasing c 2 s decreases the amplitude of fluc- 
tuations in the matter relative to the potential fluctu- 
ations. On the other hand, as we shall see in the next 
section, decreasing c 2 s also eliminates a source of CMB 



4.2. Viscous Scale and Radiation 



As discussed in 
cous scale through c: 



3.3, the effect of changing the vis- 
has little effect on the spectrum 



of matter fluctuations as long as \c 2 s \ > \c 2 is \. We show 
an example in Fig. |^a where the neutrinos in sCDM have 
been replaced with CDM as in Fig. ||. In the large-scale 
structure regime, the difference between these models 
and sCDM is small and c^ is = 1/3 in particular provides 
an excellent approximation to free-streaming neutrinos. 
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4.3. Time-Dependent Stresses and MDM 

In general, the stress parameters (w g , c 2 ff , c 2 is ) may 
all be time-dependent. An interesting concrete example 
of such a model is provided by the mixed dark matter 
(MDM) scenario where a component of HDM (e.g. mas- 
sive neutrinos) is added to the CDM. Here the equation 
of state goes from Whdm = 1/3 to as the neutrinos 
become non-relativistic. Fitting to the numerical inte- 
gration of the distribution gives 



1 



(o/Onr) 



2p1 



-l/p 



(21) 



with p = 0.872 and a nl = 6.32 x 1CT 6 /^ 2 - We can 
model its behavior as a GDM-component with w g — 
Whdm, = c 2 given by equation (||), and c 2 is = w hdm . 



In Fig. |5|b, we show that this model accurately repro- 
duces the fea tures of the MDM mo del as calculated by 
CMBFAST ( |5el,jak fe Zaldarriaga| 1996). The novelty 
of this type of model is that the ratio of the clustering 
scale to the horizon scale varies, in this case shrinking 
with time. This can have the effect of smoothing out 
the clustering feature in the matter power spectrum (c.f. 
Figs. | and |b). 

An exotic example of this type of model is the self- 
interacting dark matter candidate proposed by |Carlson] 
Ma Aacck fc Hall] (1992), where the equation of state 
passes a logarithmically decaying regime 



"'a 



31n(a/a) ' 



(22) 

2 . Here a is a constant. 
Other examples include a scalar field that rolls from 
a potential dominated regime with w g = — 1 to a ki- 
netic energy domina ted regime with w a = 1 (dec aying- 



between the radiation and matter limits of 
and w g en a 



A scenarios, see e.g. Coble, Dodclson fc Fricmar] 1997). 
The scalar field may also go from a rolling to rapidly- 
oscillating regime where subhorizon clustering can take 
place. Here one must redefine c 2 ff to be time- variable as 
well. This may occur in certain two-field models where 
the mass term can be time-dependent (see Appendix). 

4.4. GDM-Only Models 

The freedom to set the clustering scale well below 
the horizon (c 2 ff <C 1) raises the possibility that there is 
only a single component of dark matter with w g < 0, i.e. 
CDM is absent. Conventional scalar field models (e.g. 



Caldwell, Dave & Steinhardt 1998) do not allow this 



the small amplitude they possessed at horizon crossing 
(but see Appendix). 

The lack of a CDM component allows the appearance 
of interesting phenomena in the matter power spectrum. 
The main effect is that the shape parameter of CDM 
is rescaled for a given T = Sigh, because the relevant 
scale is the horizon at GDM-radiation equality given by 
equation (Il9|), i.e. 



r o ff 

r 



-3Wg/(l-3Wg) 



(23) 



possibility since perturbations could never grow beyond 



where 0,. = 4.17 x 10~ 5 /i~ 2 with the usual thermal 
history. Note that w g appear in the exponent for the 
shape paramter T c ff and hence even a mild departure 
from zero yields dramatic effects. In Fig. ^a, we show 
an example with w g = —0.04 GDM replacing the CDM 
in a model with Q tot = 1, h = 0.65, Q b h 2 = 0.0125, 
n = 1 (solid lines). This model has r e ff = 0.24 and 
closely resembles a CDM model with V = Q, m h = 0.24 
(short-dashed lines). Not only does lowering w g from 
zero help the problems of the normalization and shape 
of the CDM power spectrum, it also raises the age. 

Of course the small change from w g = to —0.04 
in the example above, only increases the age from to 
a negligible amount from 10 to 10.4 Gyrs (h = 0.65). 
If we push w g to —0.1, the age is 11.1 Gyrs but r c ff 
is too small to account for large scale structure with a 
scale-invariant n — 1 spectrum. This sort of model can 
be made viable by blue tilting the initial spectrum. In 
Fig. ^]b, we show an example with f2 t ot = 1, h = 0.65, 
fifth 2 = 0.025 and n = 1.7 (with a COBE normalization 
to a model with reionization at z = 65 in order that the 
large tilt be consistent with degree scale anisotropies) . 
It is interesting to note that the GDM power is so sup- 
pressed that even with this big bang nucleosynthesis 
baryon content, acoustic oscillations are visible in a 
critical density model. In fact, there is an interest- 
ing feature in this model at the 100/i -1 Mpc scale of 
k ~ 0.05 — 0.07ft-Mpc _1 . These models thus escape the 
constraints presented in Eiscnstcin et al.| (1997) and 
may help to explain the observed lOO/i -1 Mpc excess 
should it persist. More generally, the replacement of 
CDM with w g < GDM gives one the freedom to in- 
crease the prominence of the acoustic oscillations in the 
matter power spectrum. 

5. CMB Anisotropies 

The presence of GDM affects the CMB anisotropies 
both by its influence on the background expansion 
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Fig. 6.— GDM-only Models, (a) w g 



-0.04. Here 



an h — 0.65 model with otherwise sCDM parameters is 
shown (solid line erg = 0.60, 050/08 = 0.17) in compar- 
ison with sCDM (long-dashed line) and a CDM model 
with r = VL m h — 0.24 (short-dashed line) with a nor- 
malization artificially set to match the GDM model, (b) 
w g = —0.1. Acoustic features appear even though the 
model is a critical universe (f2tot = 1, h = 0.65) with 
big bang nucleosynthesis baryons (f^/i 2 = 0.025). The 
suppression of GDM power is counteracted by a strong 
blue tilt n = 1.7 (a s = 0.40, a 50 /(r s = 0.23). 



5.1) and on the gravitational potential perturbations 



5.2). If w g < and CDM is also present, small- angle 



anisotropics are affected only by the former since the 
perturbations that generated them entered the horizon 
well when the GDM effects are negligible. We discuss 
the case where it is absent below in § |5.3| . Large-angle 
anisotropics depend mainly on gravitational potential 
variations (the ISW effect), since primary anisotropies 
have few features that can be shifted by a change in the 
background geometry. 

5.1. Acoustic Peaks 

The acoustic peaks in the CMB depend on the pho- 
ton to baryon ratio, the expansion rate at last scat- 
tering, and the gravitational potential, all at last scat- 
tering, as well as the angular diameter distance to last 
scattering (|Hu fe White] 1996). Provided that GDM 
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Fig. 7. — Angular diameter distance, the acoustic peaks 
and the late ISW effect. Small- angle anisotropies de- 
pend on the equation of state w g through the angular 
diameter distance cLa ■ Once the angular scale is rescaled 
to the fiducial d,A of a w g = — 1 , A model, the curves are 
indistinguishable if normalized to small scales. Large- 
angle contributions arise from the ISW effect and are 
maximal in these scalar field c 2 = 1 examples. Here 
fitot = 1, O, = 0.65, h = 0.7 and tt b h 2 = 0.0125. 



contributes negligibly to the density at last scattering, 
it can only alter the peaks through the last effect. Here 
features shift in scale with the angular diameter dis- 
tance, which in a flat universe is 



d A = m- V* 



(24) 



Her e 77^, _ rj(a^) is the conformal time at last scattering 
(see |Hu fc Whitej 1997a eq. [22]-[23]). We display this 
effect in Fig. [7| where the angular scale of 4 models with 
— 1 < w g < — 1/6 are rescaled via the d A of a fiducial 
Wg = — 1 model. As w g decreases, (1a increases such 
that features are shifted to smaller angular scales. 

5.2. ISW Effect 

A more complicated effect arises from the decay in 
the gravitational potential induced by the GDM. This 
is called the late ISW effect and produces a contribution 
to Ce as 



a 



(ISW) 



dk, 



d v (^ - <$>)j e (k(r)o - ??)) 



(25) 
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Fig. 8. — Sound speed effects. Decreasing Cg ff , for a 
fixed background cosmology with w g = —1/6 and the 
same cosmological parameters as Fig. [?], decreases the 
ISW effect at large angles such that COBE normalized 
models have lower acoustic peaks. 



which appears at large angles. On small scales, the pho- 
tons can traverse many wavelengths of the fluctuation 
during the decay time and thus destroy the coherence 
of the gravitational redshifts. The cancellation is ex- 
pressed by the integral over the oscillatory Bessel func- 
tion. 

The effect is minimized if c c fr <C 1 . A concrete exam- 
ple of this is the MDM scenario where the hot compo- 
nent has w g = and an effective sound horizon today 
well below the particle horizon. Hence there is essen- 
tially no ISW contribution in this model (see Fig. ||). 
More generally, there will be a small effect if w g ^ 
due to the mild potential variation from the change in 
the equation of state (see eq. ) . 

The effect is maximized for c c g = 1 , as is the case for 
scalar field models shown in Fig. |7|. Here the potential 
decays due to pressure support of the GDM-fluctuations 
during GDM-domination. The crucial aspect is that 
the decay occurs as soon as the perturbation crosses 
the horizon so that the photons have not had sufficient 
time to cross the perturbation. The effect thus mono- 
tonically decreases as Cg ff decreases to zero. In Fig. ||, 
we display this trend in Cg. Here, we properly normal- 
ize the spectrum to the COBE detection, which corre- 
sponds roughly to £ = 10. Thus as Cg ff increases, the 
height of the acoustic peaks decreases relative to the 
ISW-boosted large-scale anisotropy. This also has the 



consequence of decreasing the normalization of the mat- 
ter power spectrum as we have seen in Fig. ^. 
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Fig. 9. — (a) Modeling radiation. Shown here are 
the anisotropies for the model where the neutrinos 
of sCDM are replaced with GDM w g = 1/3 as in 
Fig. ||a. Changing the viscosity parameter alters the 
anisotropies with c^ is = 1/3 best approximating the 
neutrinos, (b) Modeling HDM (same parameters as 
Fig. ||b). The GDM model accurately reproduces the 
features of MDM which are themselves only slightly dif- 
ferent from sCDM. 



5.3. Special Cases 

It is worthwhile to consider a few special cases to fur- 
ther illustrate the range of phenomena and show how 
more conventional candidates are recovered. As we have 



seen in §3^, radiation can be modeled through the vis- 
cous parameter c^ is . Changing c^ is has a greater effect 
on the CMB than on the matter power spectrum since 
it enters directly into the evolution of the gravitational 
potentials (see eq. fL3f| ). In Fig. [|a, we show the model 
of Figs. H and ||a, where the neutrinos in sCDM are re- 



placed by GDM of w„ 



1/3. The model with 



1 /3 yields an excellent approximation to the neu- 



trinos, whereas that with 
the sCDM results (see also 



differs by ■ 
~1995). 



20% from 



Hu et al 



Likewis e th e GDM model for a hot component pre- 
sented in §4.3 accurately reproduces the anisotropies of 
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an MDM model (see Fig. |^b). Note that this model has 
essentially no late ISW effect and differs from sCDM 
(long-dashed lines) only through the small change in 
the expansion rate and gravitational potentials due to 
the presence of the hot component at last scattering. 
MDM models have larger small scale anisotropies due 
potential decay from the radiation pressure of the hot 
component (Ma & Bertschinger 1995; Dodclson, Gates] 
& Stebbins| 1996). These results in conjunction with 
the analogous ones for the matter power spectrum pre- 
sented in Fig. H demonstrate that the GDM ansatz also 
allows one to model the full range of leading particle 
dark matter candidates. 



ties which can vary in time but not in space. These 
are the equation of state w g , the effective sound speed 
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Fig. 10.— A GDM-only Model. CDM is replaced with 
GDM of w g = -0.04 ( C 2 ff = c 2 vis = 0) as in the model of 
Fig. [6]a. Compared with the sCDM model, the acous- 
tic peaks are of higher amplitude and larger angle as 
discussed in the text. 

Finally, we show the CMB anisotropies for the sin- 
gle dark matter component model of Fig. ^|a. Mod- 
els of this type tend to have enhanced degree scale 
anisotropies. Because radiation domination is extended, 
there is more decay in the gravitational potentials due 
to radiation pressure support. This leads to early ISW 
contributions around GDM-radiation equality. Further- 
more, the sound horizon at last scattering increases due 
to the extended period of radiation domination. This 
moves the acoustic peaks to slightly larger angles. 

6. Conclusions 

We have introduced a parameterization of a gener- 
alized dark matter component based on three quanti- 



ng, and the viscosity parameter 



Combinations 



of these parameters recover all currently popular candi- 
dates for dark matter either exactly or to high accuracy, 
e.g. cold dark matter, radiation, hot dark matter, and 
scalar fields. Note that the computational costs of mod- 
eling the HDM is the same as CDM, in sharp contrast to 
a full solution to the energy-dependent Boltzmann equa- 
tion for massive neutrinos (Ma & Bertschinger 1995). 

In this general study, we have shown how the clus- 
tering scale of the GDM is not in general specified by 
its equation of state. The clustering scale appears as a 
feature in the matter power spectrum but only a weak 
enhancement of CMB anisotropies at large angles. In 
fact, the contribution is bounded by the scalar field case 



of 



cir 



= 1 for any given equation of state. There exist 



a class of models where GDM clustering has dramatic 
effects on large scale structure but no effect on the CMB 
(within cosmic variance). Conversely, altering the vis- 
cosity parameter of the dark matter affects the CMB 
anisotropies more strongly than the matter spectrum. 
An example of this behavior is provided by massless 
neutrinos. The difference between massless neutrinos 
and a perfect fluid with w g = 1/3 produces an observ- 
able difference in the CMB anisotropy spectrum. Fi- 
nally models exist where only a single species of dark 
matter is necessary in contrast to scalar field and MDM 
models where a comparable amount of CDM must also 
be present to form structure. 

While observations currently do not force one to con- 
sider the the exotic types of dark matter studied here, 
this situation may soon change with the high precision 
measurements expected from the Microwave Anisotropy 
Probe (MAP) on the CMB side and the 2 Degree Field 
(2DF) and Sloan Digital Sky Survey (SDSS) on the 
galaxy-survey side (Hu, et al. 1998). The freedom to 



alter large scale structure in relation to the CMB uncov- 
ered here may then be essential in the reconstruction of 
the cosmological model. This is especially true since the 
ambiguity introduced by the initial spectrum of fluctua- 
tions is removed once the CMB and large scale structure 
are measured at the same physical scale. Here we have 
exposed the aspects of the dark matter to which the 
CMB and large scale structure are and are not sensitive. 
These phenomenological aspects, once observationally 
determined, should aid isolation of a viable physical 
candidate for the dark matter. 
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A. Scalar Fields as GDM 

In this appendix, we demonstrate that scalar fields 
are recovered exactly as a limiting case of the GDM 
ansatz and discuss a few special cases. A minimally 
coupled scalar field ip with the Lagrangian 



-V~g[g^d^d uV + 2V{p)] 



and small perturbations tp = 



obeys 



t>o + a 2 V. v = , 



(Al) 



(A2) 



for the background field and 



-2-01 - (k 2 + a 2 V^ w )cj)i 



+ (h v - Shs)4> - 2a 2 V v h v 



(A3) 



for the perturbations. Recall that the metric perturba- 
tions h v and hg were defined for the synchronous and 
Newtonian gauges in equation (|?j). 

From the stress-energy tensor 



i 

we can associate 



1 



(V>"V;a + 2V)6» V , (A4) 



= ^a- 2 4>l + V, 
= \a- 2 4>l-V, 



(A5) 



for the background and 



{P4>+P<t>) v 4> 



b 4>i - 4>lh v ) + V v 9i i 



2V 



V <PI 



P^4, = . 



(A6) 



for the perturbations. 

The equations of motion can now be rewritten as 



= -3(^0 + P4>)~ ' 



(A7) 



for the background and 



One can rewrite the equations (|1 
GDM as 



—{p,!) +P4 > ){kv <t> + 3hs) - - 6Sp,p 

+9-(l - cDip^+p^v^/k 

-(1 + 3c|)(p +p ( f > )-v (f ,/k 
+5p^ + (p0 + p<f,)h v . (A8) 

T), and (ph for the 



(s Pg y 



[{Pg+Pg)Vg/k\ 



-{p g +Pg){kv g + 3h$) 
-3-PgFg - 3(1 + c 2 )-8p g 



-^-(pg + Pg)Vg/k + C 2 g Sp g 



-Pg T 9 - oi 1 - 3K / k2 )Pg w 9 



-{p<t> +P<p)kh v . 



(A9) 



Employing equation (||) and (|12j), we find that GDM 



with w 



-off 



-vis 



corresponds to 



a scalar field component. 

Although exact, there are a few novel aspects con- 
cerning scalar fields that are hidden in this representa- 
tion. The equation of state w g encodes the potential 
V but is not completely specified by it. The equation 
of state is also dependent on how the field rolls in the 
potential as a function of time and this is dependent on 
the expansion rate which acts a frictional term in equa- 
tion (A2). This fact allows for novel features such as 



the attractor solutions investigated by Ferreira fc Joyce 
(1997). 

Furthermore quadratic potentials V — ■^in 2 (f 2 such 
as found in axion models cause the field and so the equa- 
tion of state to rapidly oscillate. Here equation (A3) 
acts as a driven oscillator such the perturbation does not 
stabilize until (Khlopov, Malomed & Zcl'dovich 1985; 



Nambu fc Sasaki 1990) 



k 2 k, m^/Gpg, 



(A10) 



despite the fact that cjL = 1. Although the GDM de- 



scription of equation (AE) formally holds, it is impracti- 
cal to solve because of the the widely-separated expan- 
sion and oscillation time scales. One can however model 
this system by time averaging the oscillations and set- 
ting w g — and c 2 s <C 1 as is commonly done for axion 
models. 
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The field may be both rapidly oscillating and slowly 
rolling in certain two-field models (see e.g. Kodama & 
Sas |aki| 1984 for multifield modifications to the equations 
of motion). Since the field evolution depends only on 
potential derivatives, the axion-type instability exists 
on short time scales. Nevertheless on the expansion 
time scale, the rolling contributes to the background 
density and pressure. Models of this type can have a 
net w g < but can cluster well below the horizon scale. 
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